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Point interactions in a tube 

Pavel Exner 



Abstract. We discuss discuss spectral and scattering properties of a particle 
confined to a straight Dirichlet tube in ffi 3 with a family of point interactions. 



Point interactions belong to the list of problems to which Sergio Albeverio made a 
significant contribution. This topic combines a practical importance as a source of 
numerous solvable models with an aesthetic appeal as the monograph Q witnesses. 
At the same time it is far to be closed; despite the extensive and thorough character 
. of the mentioned treatise new questions still arise. 

One of them concerns point interactions in tubular regions which represent a 
natural model for a "quantum wire" with impurities. The simplest situation when 
the tube is a straight planar strip was investigated in [Q]; we refer to this paper 
for a detailed motivation and bibliography. In the present paper we present a brief 
discussion of a straight tube in R 3 with a family of point interactions. 
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1. The one-center case 



Let := 1 x M where M C R 2 is a closed compact set; we suppose that it is 
pathwise connected and dM has the segment property j^] . The free Hamiltonian is 
the corresponding Dirichlet Laplacian, Hq = — A^> with the domain Wq ,2 {Vl). It can 
, be expressed by means of the one-dimensional Laplacian and — A^. The last named 

5_i ' operator has a purely discrete spectrum; we denote by x„, v n its eigenfuctions and 

eigenvalues, respectively. For any z G C \ [^0j°°) the free resolvent is an integral 
operator with the kernel 

I e ik n (z)\x 1 -x 2 \ 

(1.1) G ($i,x 2 ;z) = (H - z) _1 (xi,x 2 ) = - V — — Xn(yi)Xn(V2) , 

1 n=0 ' tn( - Z) 



where Xj — (xj,yj) and k n (z) := ^/z — v n , which is defined and smooth except at 

x\ = x-2- It is a multivalued function of z with cuts [v n , oo), n = 0, 1, 

Suppose now that a point interaction is situated at a = (a,b) £ M°. We define 
it as in ffl], i.e. as a self-adjoint extension of the operator — A^JQj ^ \ {a}). We 
employ generalized boundary values 



(1.2) Lo(ip,a) := lim ip(x)\x— a\ , Li(ip,a) := lim__ 



L (ip,a) 
q 
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then the extension in question is specified by the boundary condition 
(1.3) L 1 (i/;,a)+<i'KaLo(i>,a) = 

for a given a € M. We shall denote it H(a,a); the case a = oo, i.e. Lq(iP,ci) = 0, 
corresponds to the free Hamiltonian Hq. 

The resolvent of H (a, a) is obtained by Krein's formula. Mimicking the argu- 
ment of |Q we get 

G (xi,a;z)G (a,x 2 ;z) 



(1.4) 



(H(a, a) — z) (xi,x 2 ) = Gq(xi, x 2 ; z) + 



a — £(a; z) 



where £(a; z) is the regularized Green's function at a, 

Ak n u 



(1.5) 



£(a; z) 



lim 



E 



\Xn{b)Y 
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The existence of the limit follows from the kernel behaviour at the singularity Q. 
However, we also need a prescription how to compute it and this differs from the 
two-dimensional case. We use semiclassical properties of the above series terms 
|| . The transverse eigenvalues behave as v n as 47r|A/|~ 1 ri for n — * oo so k n — 
2jy / 7r|M|~ 1 / 2 y / n + 0(\). On the other hand, the probability densities \xn\ 2 are 
rapidly oscillating functions. Since M supports no potential the mean value of these 
oscillations equals the constant |M| _1 ; assuming that |x n (6)| 2 oscillates around this 
value as n — > oo we can assess the divergence rate of the first series. Next we use 
the identity 

e -7«v/s 
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■ ds 



with 7 := 2y / 7r|M| and (3 1 := \\J -k\M\ and write the r.h.s. as a sum of the 
integrals over (n, n+1) obtaining 



^(a; z) = lim 



E 

n=0 



2^ n 



IXn(6)| 5 



4iru 



where K n (z) := —ik n (z) = y/i/ n — z. In combination with the preceding argument, 
it is easy to check that the summand has a uniform bound of order o(rt~ 3 / 2 ). Hence 
the limit can be interchanged with the sum and 



(1.6) 



£(a;z) 



E 

71 = 



lx»WI 2 

2n n (z) 



Vn+1 



2y/^\M\ 



Remarks 1.1. (i) We do not give here details of the oscillation argument. 
Notice that the conclusions made below can be obtained even without it, up to an 
additive renormalization of the coupling constant, since by the Weyl formula and 
the uniform boundedness of the Xn' s the difference £(a; z) — £(a; zq) is given by a 
convergent series. 

(ii) The scaling behaviour for Q a = K x M a with M a := aM, a > is more 
complicated than in the two-dimensional case. We have £(a CT ; zcr~ 2 ) = a _1 ^(a; z), 
so the singularities of the resolvent kernel are related by 

(1.7) e a {a a ,a a ) = a- 2 e{a,a) , a° 
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Proposition 1.2. The operator H(a, a) has for any a € R a single eigenvalue 
e(a,a) £ (—oo^vq). The corresponding eigenfunction is 

™ e - Kn (e)\x-a\ 

(1-8) ip(x;a,a) = \ — 5 — — Xn{y)Xn{b) ■ 

The function e(-,a) is strictly increasing and behaves as 

(1-9) e(a,3) = i*-^^ +0(a- 3 ), 

m i/ie Zzmii of weak coupling, a — » +00. Moreover, there are no eigenvalues embed- 
ded in a c (H(a, a)) = [vq, 00). 



Proof: Due to (|l.6| ), £(a; •) is strictly increasing with Ran£ = K and £(<z; 2) = 
^ lxo(fe)| 2 (^o — z)~ 1//2 + 0(1) as z — > i/q— . The non-normalized eigenfunction (1.8) is 
given by the residue term in (1.4). To check the absence of embedded eigenvalues we 
have to show that £(a; z) — a has no solutions on [vq, 00). Away of the thresholds, 
this follows from 

(1.10) Im£(a;,)= £ i^EL > . 

If |xn(fr)| 2 ^ 0, the resolvent kernel has a finite limit as z approaches v n , otherwise 
it has the same singularity as Go(xi, a; •) there, so in neither case it has a pole. | 

Mimicking the argument of [Q] we also get 

Proposition 1.3. The on-shell S-matrix at energy z — k 2 is a 2A opon x 2AT opon 
unitary matrix with elementary blocks 

/, n\ o k m I t nm r nm \ 1 at 

(1.11) b nm = (/— r , n,TO= 1, . . . , iVopen , 

V \ 'nm ^nrn J 

where N opcn := card{^„ : v n < z}, the tilded quantities are obtained by switching 
sign of the longitudinal component of a, an —a, and 

(1.12) r nm e~ lk ™ a = (t nm - S nm )e tkma = 6 Xn{b)Xm(b) ■ 

2k m a - £(a; z) 

2. Finite number of perturbations 

Denote a := {eti, . . . a at}, where Oj = (aj, bj), and a := {ai, . . . , un}, j = 1, . . . , N. 
The Hamiltonian (H(a, a) with N point interactions is defined as the self-adjoint 
extension of the operator — A^/[C^°(n \ {a}) specified by the boundary conditions 

(2.1) Li(ip, dj) + 4 7 rajL (V ) , 8j) =0, j = 1, . . . , N . 

The resolvent is again found by means of the Krein formula: 
(2.2) 

N 

(H(a,a)-zy 1 (x 1 ,x 2 ) = G (xi, x 2 ; z) + 2J Xjk(ot,a;z)Go(xi,a;z)Go(a,x 2 ;z) , 
where A(a, a; z) = A(a, a; z) _1 with 

(2.3) Ajj = otj - £,(af,z) , A jk = -G (aj,a k ; z) for j ^ k, 
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where £(dj',z) is given by ( |l.5| ), (pT^). With these prerequisites we can derive 
spectral properties of our point-interaction Hamiltonian. 

Theorem 2.1. (a) The spectrum of H(a,a) consists for any a £ M. N of the 
absolutely continuous part [vq,oo) and eigenvalues t\ < 62 < • • • < e rn < vq with 
1 < m < N , given by the condition 

(2.4) detA(a,a,z) = 0. 

The corresponding eigenf unctions are tp(x) — ^2j = idjGo(x,aj] z), where d € l w 

solves X)m=i h-{ z )jmdm = 0. The ground-state eigenfunction is positive. 

(b) z > vo cannot be an eigenvalue corresponding to an eigenvector from the sub- 
space @{ n . v <Z }L 2 (M.) <g> {Xn}- On the other hand, H(a,a) can have embedded 
eigenvalues if the family {57, a, a} has a suitable symmetry. 

(c) In the weak coupling limit, \A\ := mhii<j<jv otj — > oo, there is a single eigen- 
value which behaves as 

( f^N 



(2.5) e(a,a) = v 



2 EiLi a jXo(bj 



0(\A\ 



-3\ 



Proof: (a) A finite-rank perturbation in the resolvent preserves <J a c(Ho) = 
[uq, oo). The discrete spectrum is determined by poles of the resolvent coming from 
the coefficients Xjk in (2.2). This yields (2.4); the eigenfunctions are obtained in 
the same way as in [[I], |4[. The next question concerns the existence of solutions to 
( [2.4| ). If z — ► — oo the matrix can be written as £(a; z)A(a, a, z), where A — > —I, 
hence all eigenvalues of A(a, a, z) tend to +oo. On the other hand, for z — > vq— 
we have 

A(a,a,z) = - jL M 1 + 0(1), 

2^v - z 

where M\ := (xo(^j)xo(^m))^ m= i • This matrix has, in particular, an eigenvector 
(Xo(^i)) • • • i Xo(&iv)) corresponding to the positive eigenvalue Y^jLi Xo(bj) 2 , and 
therefore at least one of the eigenvalues of A(a, a, z) tends to — oo as z — > i/q— . 
Using the continuity we see that there is an eigenvalue which crosses zero, i.e. 
H(a, a) has at least one eigenvalue. By a straightforward differentiation we find 



— A(z) jm = - — /2 (1 + \a,j - a m yu n - z) Xo(bj)xo(b m ) ■ 

az n=Q ^{yn z ) 

The matrix function A(-) is monotonous if for any c S the quantity 4-{c, A(z)c) 
has a definite sign (is non-positive in our case). This is true provided the function 
/ : f(x) = e~ K > x >(l + n\x\) is of positive type for any n > 0, which follows from 
the identity 

O k 3 f „ipx 

(1 + K\x\)e-"W = — / , * ,„ dp 



<• JK {P 2 i - J 

and Bochner's theorem ^ Sec. IX. 2]. In fact, since the measure in the last integral is 
pointwise positive, 4z A(z) is even strictly positive; it means that all the eigenvalues 
of A(a, a; z) are decreasing functions of z and H(a, a) has at most N eigenvalues. 

To check that the ground state is non-degenerate, we have to demonstrate that 
the lowest eigenvalue of A(z) is simple for any z <G (—00,^0), which is equivalent 
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to the claim that the matrix semigroup { e -* A ( z ) : t > } is positivity preserving 
|5|, Sec. XIII. 12]. The last property is ensured if all the non-diagonal elements of 



A(z) are negative; we have A(z)j m = —Ga(aj,a m ;z) by (2.3) so the desired result 
follows from the positivity of the free-resolvent kernel. The coefficients may be 
therefore chosen of the same sign for the ground state; in fact, as strictly positive 
because dj — would mean that the cigcnfunction is smooth at x — aj so the 
corresponding interaction is absent, a JO = 00. 

(b) Suppose now that Hip — zip for some z > Vq. We adapt again the argument 
from pi Sec. II. 1] and pick an arbitrary z' € p(H); then there is a vector tpQ S D(H ) 
which allows us to write 

N 

(2.6) (p = i} +J2 d i G o{->af,z'). 

i=i 

Furthermore, we expand tpo as a series, ipo(x) — S^Lo 9n(%)Xn(y) with the coef- 
ficients g n £ L 2 (IR). Using the identity (Ho — z)ip Q = (z — z') J2f=i djGo(-, a -; z') 
and the fact that {xn} is an orthonormal basis in L 2 (M), we obtain a system of 
equations; by the Fourier-Plancherel operator it is transformed into 

z — z' - e~ ipaj 

(2-7) (p 2 - z + v n )g n (p) = J2 d iXn(bj) - - - — . 

i=i " 
If g n G L 2 the same has to be true for g n ; this is impossible if z > v n and the 



r.h.s. of (2.7) is nonzero at ±p n , where p n := \Jv n — z, since g^ would have then a 
non-integrable singularity. If N > 1, it might happen that the r.h.s. of fl2.7| ) is not 
zero identically. However, if the a,j are mutually different, Y^,j=i d jXn(bj) e Twaj — 
implies dj — by linear independence. On the other hand, if some of them coincide 
we find Y) j dj\ n (bj) — where the index runs through the values with the same 
longitudinal coordinate ay, and therefore g n — again. 

The condition v n < z in the above argument is crucial; the operator H(a,a) 
can have embedded eigenvalues with eigenfunctions in the orthogonal complement 
of the mentioned subspace if TV > 1. Examples can be constructed as in [Q (or 
other similar systems - cf.[^]) using M with a symmetry: one has to choose a family 
of weak enough point interactions with the same symmetry. 

(c) We denote A := diag(«i, . . . ,c*/v), and use the decomposition 



A(z) = (A-T(z)) \l-(A-f(z)) 



Mi 



2^0 - z 



where T(z) is a remainder independent of a , whose norm is bounded as z-> vq— , 
and Mi is the matrix defined above. The first factor is regular for |A| large enough. 
Since Mi is rank one we have to solve the equation 

N N 

v E xo&te ( 7 - r^)^ 1 ) xo(h) - E xo(b 3 f = 
j.k=i 3 , 1 



with 77 := 2^uq — z; then (2.Z) follows by the implicit-function theorem. | 
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Remarks 2.2. (i) We get also the weak- coupling asymptotics for the eigen- 
function: 



(EjLixo(&i) 2 ) i=i 



+ V Xr»(&j) V / Xn(&j')Xo(&j) • 

n=l j = l V"» - ^ 

The leading term is a product of Xi(y) with a linear combination of the eigenfunc- 
tions of one-dimensional point interactions placed at a 3 ■, j = 1, . . . , N . 
(ii) The scattering problem can be treated as in the one-center case. Existence 
and completeness of wave operators follow from the Kato-Birman theory 0] . The 
reflection and transmission amplitudes from the n-th to the m-th channel are 

rw(z) = - J2 (Hz)-i) jk Xmi l j) ff k) e *(W^> , 

(2.8) t nm {z) = 6 nm + i V (A(z)- 1 ) jfc Xro( y^ (6fc) e -^- k ^ 
and the unitarity condition now reads 

{m: Q<u m <z] 

(2.9) ^ ^ &m {inm^ sm ^nm^slri) — 5 
{m: 0<u m <z} 

because S nm is given again by (1.11), where the tilded quantities are obtained by 
mirror transformation, dj — > — Oj . 

3. The periodic case 

In the infinite-center case we restrict ourselves to the periodic situation, i.e. we 
suppose that the set {a, a} per = {[otj, Sj] : j = 1, 2, . . . } in Q is countably infinite 
and has a periodic pattern with a period £ > and N perturbations in each cell, 
which we denote again as {a, a}. Following the Floquet-Bloch decomposition, we 
find the unitary operator U : L 2 (Q) -> L 2 (B, (£/2n)d9; L 2 (Cl)), where 



x M 



(3.1) ft :=[(),£) xM, B:=[-- V j) 

the ^-projections of these sets are the Wigner-Seitz cell of the underlying one- 
dimensional lattice and the corresponding Brillouin zone, respectively. By means 
of U , the operator H (a, a) is unitarily equivalent to 

£ r® 

(3.2) U H{{a,c[) per )U- 1 = — H(a,a;0)d0, 

ZlT J\ee\<ir 

where H(a,a;9) is the point-interaction Hamiltonian on L 2 (fl), i.e. the Laplacian 
satisfying (2.1) at the points dj, Dirichlct b.c. for x E [0,£) , y £ dM, and 

(3.3) W-,ff) = e m 4'(0+, y) , ^(l- y) = e m |^(0+, y) 
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for y g M. The resolvent of H(a, o; 9) can be derived by modifying the argument 
of |Q Sec. 5]. The "free" eigenvalues 

/ 2-nm, \ 2 

(3.4) e mn (9) := ( + 9j +v n , m g Z , n = 0, 2, . . . , 

correspond to the eigenfunctions rf m ® \n, where \ n are as above and rf m {x) :— 
£-1/2 e i(2mn+0e)x/e ^ m £ z. Moreover, the free resolvent kernel is in analogy with 
JiJ obtained by a partial summation of the appropriate double series and equals 

_ ^ S mh((£-\x 1 -x 2 \)\/i'n- z) + e 2lT > ee smhdxx- x 2 \\/v n - z) 
G ( Xl ,xr,e;z) - ^ cosh^n-*) - cos(fll) 

^3 w Xn(yi)Xn(V2) 



2^/v n -z 

where r\ := sgn(xi— x 2 ). The full kernel is then expressed by a formula analogous 
to (U) with 

(3.6) X(a,a, 9; z) — A(a, a, 9; z)^ 1 , 

where Aj r = —Go(3j,a r ) for j ^ r, while the diagonal elements are given by 



1 ■ ' " j 2 ^ ^cosh(^V^Fi)-cos^ ~z ^\M\ J ' 

We may also write the last formula as Ajj(a, a, 9; z) — oij — 9; z), where the 
function £ is for z£l more explicitly given by 



= 1 -sp / sin(£y z-v n ) Xnjbj) 2 y/n - y/n+l 
,Z ' ~ 2 {n .^ <2} ^cos(VF^)-cos^ V^ 3 ^ v^P*l 



1 / sinh(£Vz/»-z) Xnjbj) 2 y/n-y/n+1 

2 {n .7^ >2} ycosh^v^Fi) -cos6>£ y/U^ x/tW\ 

It is defined everywhere except at 

(3.9) 8(a, 9) := { e m „(0) g £(0) : x„(&,) ^ } , 



where E{9) is the eigenvalue set (3.4). Of course, the r.h.s. of ( |3.9|) makes no sense 
if z = u n and Xn(bj) 7^ 0, but applying general results on self-adjoint extensions 
to the one-center case, we can establish a posteriori that £ can be defined there by 
continuity. Moreover, we find that the function £(Sj, 9; •) is monotonously increasing 
between any pair of neighboring singularities. Then we have the following result. 

Proposition 3.1. For a given a g R N the operator H(a,a;9) has N eigen- 
values e mn j{a, a; 6), j = 1, . . . , N, in any gap of the set (3.1) determined by 

(3.10) det A(a, a, 0; z) = 0. 

The corresponding eigenfunctions are ip(x) — Y]f—i djGo(x, Sj, 9; z) where the dj 's 



are determined by A(a, a, 9; z) as in Theorem \2. 1\ 

The spectrum of the original Hamiltonian H({a, a} per ) consist then of bands, 
(3.11) a(H({a,a} per )) = |J {e mnj [a,a;6) : 9 g B} . 
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One is interested, of course, in its absolute continuity and existence of gaps. We 
restrict ourselves to the simplest nontrivial situation. 

Example 3.2. L et N = 1, i.e. let each cell contain a single point interaction. 
The condition ( [3.10 ) then simplifies to 



(3.12) £{S,6;z)=a. 

The left hand side is monotonously increasing between its singularities, i.e. the 
points of £ [a, 9). This means that for fixed a, 9 there is a sequence {e r (o:, a, 9)}^L 
arranged in the ascending order; each of them depends, in fact, only on the trans- 
verse component b of the vector a. The lowest one satisfies 

e (a, a, 9) < i> + 9 2 

and between each two neighboring points of £(a,9) there is just one of the other 
eigenvalues. It is also clear that any of e r (a, a, 9) is continuous with respect to 
the parameters and e r (-,a, 9) is increasing for fixed b and 9. Concerning the 9- 
dependence, the implicit-function theorem tells us that 

de r (a,a,9) 8£(a,9;z) ( d£(a,9; z)\ 1 



89 89 \ dz 

(e T ,e) 

whenever the denominator is nonzero. Away of the thresholds, z = v n , and the 
points of £(9), a straightforward differentiation shows that £(a, •; •) is analytic in 
both variables. By ( |3.8| ) the numerator is not identically zero; hence the derivative 
8e r (a, a, 9) /89 may be zero at some points but never in an interval and the spectrum 
of H({a,a} per ) is absolutely continuous |5|, Sec. XIII. 16]. 

Let us turn now to the question about the number of gaps. Below z = vq 
the spectrum may be estimated by means of extrema of the function £ which yield 
^-independent bounds: we have £(a, 9;z) < £+(a, z) where 



£+(a, z) := max £(a,9:z) 
\ee\<7T 



2 ^ 

n=0 



XoW 2 

\fVn - Z 



coth 



fn- y/n+l 



and a similar formula for the minimum, £_ (a, z), with coth replaced by tanh. Both 
functions are continuously increasing and tend to — oo as z — * — oo. On the other 
hand, ^+(a, •) diverges as z — > isq— while the lower bound ^-(a, •) has a finite limit. 
This shows, in particular, that the spectral condition ( 3.1 2| ) has no solution for any 
8 in a left neighborhood of vq provided 



< f_(a,Mo-) 



(3.13) 



1 OO 



tanh -\/v n — v + 



Vn+T 



VMM] 



in other words, that a gap exists. The condition (3.13) is satisfied for a strong 
enough coupling if the point-interaction spacing is kept fixed. On the other hand, 
inspecting the right-hand-side we see that the gap exists also for any fixed a and the 
spacing I large enough. In this respect the spectrum is similar to that of a straight 
polymer in R 3 described in |ll Sec. III. 1]. However, for our "coated polymer" a much 
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stronger result is valid: we shall show that under a suitable choice of parameters it 
can have any finite number of gaps. 

To this end, we consider z 6 (^o + £, V\ — e) for a fixed e > and £ ^ y\M |, 
and rewrite the right hand side of the relation (3.8) as 



£(a,9;z) = £ (a,6;z) +rj(a,i 

where 



£ (a,6;z) — 



cos(€\/z — isq) — cos 91 2^/z — vq 

and rj(a, 9; z) is the rest. The latter is monotonously increasing with respect to z 
and its derivative is bounded everywhere below the second threshold, in particular, 
in the chosen interval of energies. Moreover, rj{a, 9; z) is bounded from above by 



rj + (a,z) := + - > — i= coth -\/v n 



and the corresponding minimum, r)-(a,z), is obtained when coth in the last ex- 
pression is replaced by tanh. These estimates shrink as i becomes large: using the 
inequality cothu — tanhit < 5e~ 2u for 2u > 1, we find 

/-» \ /-. n . 5 Xi(&) 2 -i^=i . 5 ^ Xn{b) 2 -t^-~z 
r 1+ (a,z)-r 1 ^a,z)<-- / ==e + _ ^ __ e ^" 

The series can be estimated by an integral, which yields for a fixed z € (vq+e, v\ — e) 
the behavior 

(3.14) V+(a,z) -r)_(a,z) = O (x/\M\ • 

On the other hand, the function gg (u) := sinu(cosw— cosOl)^ 1 is increasing between 
any two zeros of its denominator. In the intervals, where it is positive, it can 
be estimated from below by the appropriate branch of tan (§ + 7rm); when it is 
negative, we have a similar estimate from above with tan replaced by — cot. Hence 
independently of 9 we have either 

Co {a,6;z) > tan - <^ -y/z~v 

or 

€o{a,9;z) < - cot - < -Vz-u 

2y/z-v \2 [it 

where {•} denotes the fractional part. Putting the estimates together we see that 
the oscillating part dominates, so for sufficiently large \a\ there are gaps having 

+ ( IL Y L ) 2 as one endpoint provided it belongs to (v a + e, V\ — e). In addition, 
tan u + cot u > 2 which means that the gap between the lower and the upper bound 
to £(a, 9; z) never closes within the chosen interval if Ij \J\M\ is large enough; we 
infer that for any a € R the operator ({a, a} pe r) can have an arbitrary finite 
number of gaps in its spectrum provided the spacing of the point-interaction array 
is large enough. 
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Conclusions of the example can be extended to a finite number of point per- 
turbations per cell. In a similar way one can treat a toroidal tube supporting point 
interactions and threaded by a magnetic flux, etc. 

A more difficult question concerns the finiteness of the number of open gaps. 
Recall that the gap number depends strongly on the dimension of a periodic system: 
for one-dimensional systems it is generically infinite J5], Sec. XIII. 16], [jj], Sec. III. 2], 
while for higher dimensions it is finite by the Bethe-Sommerfeld conjecture. The 
latter is known to be true, in particular, for periodic potentials or lattices of point 
interactions in R 3 - see [| . The tube boundary makes things more complicated, 
but one still expects that at high energies gaps will close due to overlapping of 
contributions from different transverse modes. Notice in this connection that in 
the example we have been looking for gaps in the energy interval where transport 
is possible in the lowest transverse mode only. Nevertheless, it is not easy to 
demonstrate that no gaps remain open above a certain energy. 

In a similar way, the "mixed dimensionality" of waveguide systems inspires 
other questions such as existence of a mobility edge in tubes with random point 
interactions, etc. Generally speaking, solvable models whose genealogy can be 
traced back to the treatise Jl| will represent for long a useful laboratory for the 
spectral theory. 
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